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A real space Renormalization Group approach is presented for a non-mean field spin-glass. This 
approach has been conceived in the effort to develop an alternative method to the Renormalization 
Group approaches based on the replica method. Indeed, non-perturbative effects in the latter are 
quite generally out of control, in such a way that these approaches are non-predictive. On the 
contrary, we show that the real space method developed in this work yields precise predictions for 
the critical behavior and exponents of the model. 
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I. INTRODUCTION 

Spin-glasses, structural glasses, and the physical de- 
scription of their critical properties have interested sta- 
tistical physicists for several decades. The mean-field 
theory of these models [6, 14] provides a physically and 
mathematically rich picture of their physics and of their 
critical behavior. Notwithstanding the great success of 
such mean-field theories, real spin-glasses are non-mean 
field systems, because they have short-range interactions. 
It follows that these systems cannot be described by 
mean-field models. Indeed, the generalization of the 
above mean-field theories to the non-mean field case is an 
extremely difficult task that has still not been achieved, 
so that the development of a predictive and consistent 
theory of glassy phenomena for real systems is still one 
of the most hotly debated and challenging problems in 
this domain [9]. 

This task is difficult to achieve because the perturba- 
tive field-theoretical techniques [16, 17] yielding the Ising 
model critical exponents with striking agreement with ex- 
perimental data do not apply to locally-interacting glassy 
systems. Indeed, a considerable difficulty in the set-up of 
a loop-expansion for a spin-glass with local interactions, 
is that the mean-field saddle-point has a very compli- 
cated structure ]7, 14], and non-perturbative effects are 
not under control, in such a way that the properties of 
the large scale behavior of these systems are still far from 
being clarified. 

The physical properties of the paramagnetic- 
ferromagnetic transition emerge clearly in ferromagnetic 
systems in the original work of Wilson [16], where 
one can write a simple Renormalization Group (RG) 
transformation, describing a flow under length-scale 
re-parametrizations. Later on, people found out that 
these RG equations are exact in non-mean field models 
with power-law ferromagnetic interactions, built on hi- 
erarchical lattices as Dyson Hierarchical Model (DHM) 
]8]. Indeed, in these models one can write an exact RG 
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transformation for the probability distribution of the 
magnetization of the system, in such a way that all the 
relevant physical information on criticality and all the 
fundamental RG concepts are encoded into this equa- 
tion, whose solution can be explicitly built up with the 
e-expansion technique [4, 5]. Accordingly, to investigate 
the RG properties of non-mean field spin-glasses it is 
natural to consider spin-glass models built on hierarchi- 
cal lattices. This study has been done heretofore only 
for some particular models. On the one hand, models 
with local interactions on hierarchical lattices built on 
diamond plaquettes [1[, have been widely studied in 
their spin-glass version, and have been shown to lead to 
weakly-frustrated systems even in their mean-field limit 
[12[, and so are not a good representative of a realistic 
strongly-frustrated spin-glass. On the other hand, a 
RG analysis of a different kind of random models on 
Dyson hierarchical lattice, and of their physical and 
non-physical infrared (IR) fixed points, has been done 
heretofore [15[. Unfortunately, also in these models spins 
belonging to the same hierarchical block interact each 
other with the same [15[ random coupling J, in such 
a way that frustration turns out to be relatively weak 
and they are not a good representative of a realistic 
strongly-frustrated system. 

The study of such non-mean field strongly-frustrated 
spin-glasses is difficult, also because it is hard to identify 
the correct order parameter and write the resulting RG 
equations for a function or functional of it without relying 
on the replica method, which is generally able to make 
predictions for the critical exponents only in the mean- 
field case [7[. 

In this letter we present a real space RG method for 
a non-mean field strongly-frustrated spin-glass on hierar- 
chical lattice, the Hierarchical Edwards- Anderson model 
(HEA) [11[, that does not rely on the replica method. 
Even if this method does not identify the order parame- 
ter of the system, it is interesting from the methodolog- 
ical point of view, because it yields a way to implement 
Kadanoff's [13] block-spin decimation rule in a strongly- 
frustrated system, and to write the resulting RG equa- 
tions. In this way, precise predictions on the critical ex- 
ponents are obtained. 
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The HEA is defined as a system of 2*^+^ spins 
Si, . . . , S'2fc+i, Si ~ ±1, with an energy function defined 
recursively by coupling two systems, say system 1 and 
system 2, of 2*^ spins 

^^fe+l [Si, ■ ■ ■ , S2k + l] = Hi [Si, . . . , S2k] + (1) 

+Hl[S2k+i, . . . , S2k+i] - 2^"'^''^^^^^ ^ JijSiSj, 

i=l j=2''+l 

where are random couplings distributed according to 
a Gaussian law with zero mean and unit variance, and 
Hq [S] =0. CT is a parameter tuning the decay of the inter- 
action strength between spins with distance. It turns out 
that for CT < 1/2 the thermodynamic limit is ill-defined, 
because the interaction energy grows with k faster than 
the volume 2^, while for ct > 1 the interaction energy 
goes to zero for large k, and no finite-temperature phase 
transition occurs. Accordingly, in the following we will 
take 1/2<CT<1. In this interval, the model is a non- 
mean field one, and the mean-field limit is recovered for 
CT ^ 1/2 [11]. 

The critical properties of the HEA have been stud- 
ied heretofore within the replica formalism [3], showing 
that the system has a classical behavior in the region 
1/2 < CT < 2/3, where the mean-field approximation is 
correct, while non-mean field effects are important for 
2/3 < CT < 1. This analysis makes a prediction for the 
critical exponents only in the classical region, since in the 
non-classical region the first few orders of the ct — 2/3 = e- 
expansion have a non-convergent behavior, and higher 
orders are not known. 

Before exposing the real space approach for the HEA, 
let us illustrate it in the case where the couplings Jij in 
eq. (1) are ferromagnetic, i. e. in for the well-know DHM 
[8], in order to test of the consistency of our method. 



II. THE REAL SPACE APPROACH FOR 
DYSON HIERARCHICAL MODEL 

DHM is defined [8] as a system of 2*''+^ spins 
Si, . . . , S2k+i, Si = ±1, with an energy function defined 
recursively by coupling two systems, say system 1 and 
system 2, of 2'' spins 

Hk+i [Si, ■ ■ ■ , S2k+i] = H[ [Si,. . ., S2k] + (2) 

+Hk [S2''+l, • • • ) 5*2^ + 1] + 

-j2^(--)(^+^)f-i,x:^^.l ' 



where Hq[S] = 0, the suffix F stands for 'ferromagnetic', 
and one can show that 1/2 < ap < 1, with the same 
considerations as those used to derive the constraints on 
CT for the HEA. 

The real space RG method is built up by iterating ex- 
actly the recursion equation (2) for k = ko steps. In 



this way, a DHM model with 2'=" spins >S'i, . . . , 5'2fco and 
Hamiltonian H^^[Si, . . . , 5'2fco] is obtained. We now want 
to build up a 2'^°+^-spin DHM starting from such a 2*^°- 
spin DHM, which can be done as follows. We consider 
a 2*''o^^-spin DHM, where J is replaced by another cou- 
pling J'. Such a 2'^''~^-spin DHM is defined by iterating 
ko — 1 times eq. (2) with J J' , and its Hamiltonian is 
H'^^^_i[S'i, ■ ■ ■ , S'^ko-i]- Given J, the coupling J' is cho- 
sen in such a way that the 2'°"~^-spin DHM represents 
as well as possible the 2'^"-spin DHM, as qualitatively 
depicted in fig. 1. This approximation is practically im- 
plemented by considering a physical observable O^^ iPJ) 
of the 2'="-spin DHM, and an observable 0^^^j^{(3J') of 
the 2'^''~^-spin DHM, where f3 is the inverse temperature. 
The normalized magnetization on the left half of the 2'^"- 
spin DHM is 



rriL 



2fco-l ^ 



and so for the right- half magnetization m^, where 
stands for the thermal average at fixed /3, performed with 
weight cxp(— The normalized magnetization on 
the left-half of the 2'=^o-i-spin DHM is 



2'=o-^ \ r r / 1 2*^0-^ 



2*00-2 



and so for the right-half magnetization m^. where 
E^, stands for the thermal average with weight 

exp(— /377'^^_]^). Mimicking Kadanoff's block-spin rule, 
for the 2'^°~^-spin DHM to be a good approximation of 
the 2'^''-spin DHM, we map the block of the spins in the 
left-half of the 2'^''-spin DHM, into the block of the spins 
in the left-half of the 2''«"i-spin DHM, and so for the 
right-half. To do so, we choose the observables to be 
0^,(/3J) ^ ^g[mj^mn\, 0^^_,{PJ') E^-, [m'^m'^}, and 
impose the equation 



<(/3^) = Ofo-i(/3J')- 



(3) 



For any fixed J, eq. (3) determines J' as a function of 
J, as the value of the coupling of the 2'^''^^-spin DHM 
such that this yields the best possible approximation of 
the 2'=o-spin DHM. 

Let us take two copies of the 2'^'"~^-spin DHM. We put 
these two copies into interaction according to eq. (2), 
and form a 2'^" -spin DHM. Since each of the DHMs that 
we put into interaction represents a 2'^" -spin DHM, the 
2*^" -spin DHM result of this composition effectively rep- 
resents a 2''°"'"^-spin DHM. Once again, this DHM is then 
approximated as a 2'°°^^-spin DHM with coupling, say, 
J", and so on. Setting Jq = J,Ji = J' , J2 = J",..., 
eq. (3) establishes a relation between and Jk+i, phys- 
ically representing the RG flow of the coupling Jk under 
reparametrization of the unit length 2'^ 2*^+^. 
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FIG. 1: Approximation of the real space approach for fco = 2. 
In the implementation of the real space approach to DHM, 
a 2'^-spin DHM is approximated by a 2-spin DHM. In the 
implementation of the real space approach to HEA, a 2^-spin 
HEA is approximated by a 2-spin HEA. 
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The RG eq. (3) is not exact, because it relies on the 
fact that a 2'^''-spin DHM is approximated with a 2*^""^- 
spin DHM. Nevertheless, such an approximation must 
become asymptotically exact in the large koAivciii^ where 
both 2*^" and 2'^" tend to infinity. Quite large values of 
fco can be reached by exploiting the hierarchical structure 
of the system [8], in such a way that the observables 
' ^ko-i '^^'^ calculated with a computational effort 
of 2^°. It is possible to show that for any feg the real 
space method reproduces the constraints 1/2 < CTi^ < 1. 
Indeed, for ap > 1 eq. (3) gives J' < J\/J,[3, so that the 
coupling Jfc goes to for large A:, and no phase transition 
occurs. On the contrary, for <Tp < 1/2 one has J' > 
J\/J,l3, and the model is thermodynamically unstable. 

The critical exponent vp related to the divergence 
of the correlation length [17] is easily obtained by lin- 
earizing the transformation (3 J jSJ' in the neighbor- 
hood of the critical fixed-point /3J = [3 J' = Kc [16], 
21/up ^ ^ n _ 2 depict 2^/"" computed 

with this method, together with 2^/'^^ computed by Ble- 
her [2, 4] with an independent approach, as a function 
of 1/2 < ap < 1. The latter calculation makes an ex- 
act prediction for 2^/"^^ in the region 1/2 < ap < 3/4 
where the mean-field approximation is exact, while it es- 
timates 2^/"^^ in the non-mean field region 3/4 < ap < 1 
by means of a series of successive approximations. In the 
bottom inset of fig. 2 we show how Ap fis for finite ko 
has been extrapolated to the fcp 00-limit: for every 
ap the sequence Ap us vs ko is fitted with a function 
of the form a — b ■ , and a is the resulting extrapo- 
lated value. The parameter 7 < 1 is an indicator of the 
speed of convergence with respect to k^: the larger 7 the 
slower the convergence. In the main plot of fig. 2 the 
extrapolated value is depicted, and this is in good agree- 
ment with the value given in [2]. The region where the 
disagreement between the two methods is maximum is 
ap ~ 3/4, where 2^/"^^ must be non-analytic ]5]. This 
non-analyticity cannot show up for finite ko. However, 
in the top inset of fig. 2 we show how the parameter 7 
has a maximum at ap w 3/4. This fact shows that the 
convergence slows down in the proximity of ap = 3/4, i. 
6. hat the real space method signals the appearance of a 
non-analyticity of vp at ap = 3/4, which results from the 



a fmiction of ctf for 1/2 < ap < 1. 
The black dashed line represents the upper critical dimen- 
sion ap = 3/4 discussed in [5]. The points are given by 2^''"^ 
computed with the real space method for 5 < fco < 12, and 
the gray points are 2^^'^^ extrapolated to ko — > 00 by fitting 
Af RS vs fco with a function of the form /(fco) = a — b ■ 7*°. 
The orange dashed curve and the green triangular points are 
2^'''^^ obtained in [2]. Top inset: 7 vs ap. Bottom inset: 
Afas vs fco for ap = 0.92, its fitting function /(fco) and the 
extrapolated value a. 



switchover from a mean-field to a non-mean field regime. 

It is now natural to generalize this real space approach 
to the HEA model, to compare its predictions with those 
obtained with the replica method. 



III. THE REAL SPACE APPROACH FOR THE 
HIERARCHICAL EDWARDS-ANDERSON 
MODEL 

Let us now illustrate how to apply the real space 
method to the HEA, by considering first the simplest 
case ko — 2. The reader should follow our derivation 
in close analogy with that exposed above for DHM. A 
HEA model with 2^ spins 5*1 , • • • ,6*4 and Hamiltonian 
i/2 [<S'i , . . . , 5*4] is built up exactly by means of the recur- 
sion equation (1). We set Jij = 2~'^Jy , where by defi- 
nition the couplings are independent identically- 
distributed random variables, and the probability distri- 
bution of each of them will be denoted by p{J) . Thus, we 
consider a 2-spin HEA model, whose Hamiltonian reads 
i/'i[S'J, S'2] = ~J'l2^i^2- For each realization of the cou- 
plings {>Jij}ij, we choose a function of {J^ij}ij in 
such a way that the 2-spin HEA yields the best possi- 
ble approximation of the 2^-spin HEA, as qualitatively 
depicted in fig. 1. 

In order to do so, let us consider a physical observ- 
able 02{{PJij}ii) of the 2^-spin HEA, depending on the 
6 couplings {Jij}ij and /3, and an observable Oi(/3Ji2) 
of the 2-spin HEA. Inspired by the fact that the order 
parameter in the mean- field case is the overlap [14], here 
we build up O2 and Oi as the thermal average of prod- 
ucts of spin overlaps. To build up O2 and Oi, consider 
two real replicas S^ , S^ of the spins of the 2^-spin model. 
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and two real replicas S' , S' of the spins of the 2-spin 
model. The normalized overlap between and S'^ on 
the left leaf of the 2^-spin HEA is 



Qi 



Sn So 



and so for the right-leaf overlap Qn, where stands 
for the thermal average at fixed disorder and 

p. The normalized overlap between S' and 5*' on the 
left leaf of the 2-spin HEA is Q'^ = S'^S'^, and so for 
the right-leaf overlap Q^. Following Kadanoff's deci- 
mation rule, we map the 2^-spin HEA into the 2-spin 
HEA by imposing that the spins 5*1 , ^2 correspond to 
the spin S'l , and that the spins S3 , ^4 correspond to the 
spin 5*2 . This mapping results in a correspondence be- 
tween Ql and Q^, and between Qn and Q'j^. By choos- 
ing the observables as 02({/3Jij}) = K§i s2[QlQr], 
Oi(/3Ji2) = g,2 [Q'lQ'r]i Kadanoff's decimation rule 
can be practically implemented by imposing the equality 



02({/3Jij}) = Oi(/3Ji'2): 



(4) 



where E^, stands for the thermal average at fixed dis- 
order J7]'2 and /3. For any realization of the couplings 

eq. (4) determines J7i'2 as a function of 
in such a way that the 2-spin HEA yields the best possi- 
ble approximation of the 2^-spin HEA. Accordingly, the 
distribution p{J) induces a distribution of J7i'2, that we 
will denote by v'{J\2)- The mapping between p(J') and 
p'{J^') can be shown to be given by 



(5) 



<5 ( J' - - arctanh ( \l 02[{[i J^J)) 



-5[J' + ^ arctanh ( V02({/3Jy }) 



According to the iterative construction of eq. (1), a 
new HEA is then constructed by taking two realizations 
of the 2-spin HEA. Each realization is given by throwing 
the coupling J' according to its probability distribution 
p'{J'). We put these two copies into interaction to form 
a 2^-spin HEA. Since each of the HEAs that we put into 
interaction represents a 2^-spin HEA, the 2^-spin HEA 
result of this composition effectively represents a 2'^-spin 
HEA. At the next step of the iteration, this 2^-spin HEA 
is again approximated as a 2-spin HEA with coupling, 
say, J7i'2, and the probability distribution p'\J\2) of Jyi 
is computed fromp'(j7"'), and so on. This step is repeated 
fc-times, and a system representing a 2^+''"-spin HEA is 
obtained. 

Setting po(J) = Vi{J) = P'{J). P2{J) = 

p"{J), ■ ■ ■, eq. (5) establishes a relation between Pk{J) 



and pk+i{J), physically representing the RG flow of 
the probability distribution of the coupling p^ {J) under 
reparametrization of the unit length 2^ — > 2'^"*'^. 

Eq. (5) has been solved by means of the population 
dynamics algorithm. In population dynamics, one repre- 
sents the function p{J) as a population of P 1 num- 
bers {Ji}i=i....^p, where each jTi has been drawn with 
probability distribution p{Ji). The mapping p{J^) — )■ 
p'{J') given by eq. (5) results into a mapping between 
{Ji]i and the population {Jl}i representing p'{J'). 

The structure of the fixed-points of eq. (5) has been 
thus investigated numerically, showing that there exists 
a finite value of /3 = /3c such that for /3 < /3c Pk{J) 
shrinks to a 5{J) as k is increased, while for /3 > /3c 
Pk{>J) broadens, i. e. its variance is an ever- increasing 
function of k. The physical interpretation of these two 
temperature regimes is that for /? < /3c Pki^J) flows to 
the attractive high-temperature fixed-point with J ~ 
where spins are decorrelated, while for /? > /3c it flows to 
the attractive low-temperature fixed-point with J' ~ 00 
where spins are strongly-correlated. This fact implies 
that as the temperature is lowered below Tc = l//?c 
a phase transition occurs, resulting in the appearance 
of a collective and strongly-interacting behavior of spins 
in the low-temperature phase. The existence of such a 
finite-temperature phase transition for a diluted version 
of HEA model has already been established heretofore in 
MC simulations by means of finite-size scaling techniques 
[11]. 

The population dynamics approach reproduces the fact 
that for (T < 1/2 the thermodynamic limit is ill-defined, 
which has been discussed above. Indeed, the numerics 
show that for (T — >• 1/2 /3c — > 0, in such a way that 
the variance of pk{J) , and so that of H2, is an ever- 
increasing function of fc, and the thermodynamic limit 

— > 00 is ill-defined. Unfortunately, the second con- 
straint (7 < 1 is not reproduced. This is presumably due 
to the fact that eq. (5) implements only the lowest-order 
approximation of the real space method, fco — 2, and that 
the method is exact only for large fcg- This hypothesis is 
supported by the estimate of the critical exponents that 
we will give in the following, suggesting that the closer a 
to one, the larger the values of /cq needed to have a good 
estimate of the exact result. Accordingly, for cr —>■ 1 a sig- 
nificantly better description would be obtained if larger 
values of fco were accessible, and the a < 1-limit would 
be recovered. 

The numerical implementation of eq. (5) also reveals 
the existence of a repulsive critical fixed-point with a fi- 
nite width, that we will denote by p^ {J') , which is reached 
by iterating eq. (5) with (3 = f3c. The critical exponent 
v governing the power-law divergence of the correlation 
length at /3 = /3c is determined [16] from the spectrum 
of the matrix linearizing the transformation (5) in the 
neighborhood of p*(J'). 

Before discussing the numerical results for p^{J) and 
V, let us discuss better implementations with fco > 2 of 
this method. The only new element with respect to the 
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FIG. 3: Fixed-distribution p, as a function of J for fco = 3 
and a — 0.621, 0.678, 0.737 (in black, violet and brown respec- 
tively), and for fco = 4 and cr = 0.621, 0.678, 0.737 (in red, blue 
and green respectively). For all the three values of a, the dis- 
crepancy between p*{J) in the ko — 3-approximation and 
p* ( J7) in the ko — 4-approximation is relatively small, signal- 
ing that fco = 4 is presumably large enough for the method 
to give a reasonably good estimate of the critical fixed-point. 
Inset: p.(0) vs a for ko = 3,4. A plausible picture resulting 
from the data is that, for large ko, p*(0) — for a < 2/3 and 
p*(0) > for a > 2/3. Such picture has a clear physical inter- 
pretation given in the text, and suggests a change of behavior 
at (T = 2/3, reminiscent of the switchover from a mean-field 
regime for cr < 2/3 to a non-mean field regime for a > 2/3 
predicted by the replica approach. 



ko = 2-case is the following. For /cq > 2, a 2 "-spin 
HEA is approximated as a 2'°°~^-spin HEA. The latter 
has 2'=«-i(2'=o-i - l)/2 = M' > 1 couplings {Jlj}ij. It 
turns out that even if the couplings {Jij}ij of the 2*^°- 
spin HEA are independent, {Jlj}ij are not, and are dis- 
tributed according to a joint distribution that we denote 
by Pcii'^ij}ij)- other words, correlations are intro- 
duced when iterating once the RG transformation. In 
the present treatment such correlations have been ne- 
glected by assuming that each of the {Jlj}ij behaves as 
an independent random variable distributed according to 
a distribution obtained as the average of M' distribu- 
tions, each obtained by marginalizing p'c{{'Jij}ij) with 
respect to M' — 1 couplings Jl^ . 

The real space approach has been thus implemented for 
/co = 2, 3, 4. Larger values of fco were not accessible, since 
the computational effort scales as 2^ " . All the qualita- 
tive features emerging for = 2 and discussed above, 
are preserved for /cq 3, 4. In fig. 3 we depict p^,{J) as a 
function of S for several values of a in the feg = 3, 4 ap- 
proximations. Two interesting features emerge from fig. 
3. Firstly, the discrepancy between p<,{J) in the fco = 3- 
approximation and p*{ J) in the fco = 4-approximation is 
relatively small, signaling that fco = 4 is hopefully large 
enough for the real space approach to give a reasonably 
good estimate of the critical fixed-point, at least for the 
values of a considered in fig. 3. Secondly, a plausible 
scenario resulting from the inset of fig. 3 is that, for 
large fco, p*(0) = for cr < 2/3, while p*(0) > for 



(7 > 2/3. Interestingly, the analysis of the HEA based 
on the replica approach [3, 11] predicts a sharp change of 
behavior from a mean- field regime for 1/2 < a < 2/3 to 
a non-mean field regime 2/3<tT<l. In the real space 
approach J^ij is nothing but the effective coupling be- 
tween spins Si and Sj of a 2*^" -spin HEA. At the critical 
point, Si is obtained as the coarse-graining of a group 
of 2',Z ^ 1 spins, which have been progressively deci- 
mated and reduced to a single, effective degree of free- 
dom Si , and so for Sj . In the mean- field case cr < 2/3 the 
model is mean-field, and should thus behave as a fully- 
connected one. Accordingly, the 2'-spins represented by 
Si must interact with all the other spins, and so with 
the 2'-spins representing Sj. Thus, the effective coupling 
between Si and Sj cannot vanish, i. e. p*(0) = 0. In 
the non mean- field case cr > 2/3 the system is not fully- 
connected, since the effective interaction range is finite. 
Accordingly, there is a finite probability that 2'-spins rep- 
resented by Si don't interact with the 2'-spins represent- 
ing Sj. Thus, the effective coupling between Si and Sj 
can vanish, i. e. p*(0) > 0. According to this argument, 
this change of behavior of (0) at cr = 2/3 can be seen as 
the switchover from a mean-field behavior to a non-mean 
field one, and is predicted independently and confirmed 
by the replica analysis of the HEA. 

Let us now consider the predictions on the critical ex- 
ponent v. In fig. 4 we depict 2^/"^ obtained with the 
fco = 2, 3, 4-approximation and 2^^" from the replica ap- 
proach [3] as a function of cr, both in the mean-field region 
cr < 2/3 and in the non- mean field region cr > 2/3, where 
the first two orders of the e-expansion are depicted. The 
agreement between 2^^'^ computed with the real space 
approach for fco = 2 and 2^^" computed with the replica 
approach is not satisfying. Nevertheless, for fco = 3, 4 the 
agreement in the mean- field region 1/2 < cr < 2/3 be- 
comes very good, and serves as an important test of the 
real space method. A quantitative comparison between 
2^/^^ of the real space approach and that of the replica 
approach in the non-mean field region cannot be done, 
because in the latter the e-expansion is out of control, 
i. e. the first two orders of the expansion have a non- 
convergent behavior, and higher orders are not known. 
Accordingly, the e-expansion curve depicted in fig. 4 
must not be considered as an estimate of 2^^'^ . A predic- 
tion for v in the non-classical region cr > 2/3 for a diluted 
version [11] of the HEA is given by Monte Carlo (MC) 
simulations ]10]. According to such a numerical work, 
2^^" is a decreasing function of cr in the neighborhood 
of cr = 2/3, which is in disagreement with the results of 
the real space approach, fig. 4. This discrepancy will be 
discussed in the following. 



IV. DISCUSSION AND CONCLUSIONS 

In this letter we developed a real space RG approach 
for a non-mean field spin-glass, the Hierarchical Edwards- 
Anderson model (HEA) . This approach is innovative with 
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FIG. 4: 2^''" as a function of a. The red, brown and blue 
points are 2^''" computed with the real space approach for 
fco = 2, 3, 4 respectively. The black dashed curve and the 
green dashed curve are 2^''" obtained with the replica ap- 
proach [3], and the orange dashed line is the relative upper 
critical dimension a = 2/3 [11]: the black dashed curve is the 
mean- field value of 2^^" for cr < 2/3, while the green dashed 
curve is the two-loops result obtained with the e-expansion. 

respect to the RG approaches to disordered, strongly 
frustrated systems developed heretofore that generally 
rely on the replica method [7]. Indeed, the present ap- 
proach does not make use of the replica method, which 
is generally predictive only in the mean-field case, and 
cannot handle perturbatively fiuctuations around the 
mean-field saddle-point, because these turn out to be 
out of control [7]. Through a systematic approximation 
scheme, the present approach implements Kadanoff's 
block-spin decimation rule [13] on spins. The implemen- 
tation of such a decimation rule to a disordered, strongly- 
frustrated system has not been developed heretofore be- 
cause of the intrinsic difficulties introduced by frustra- 
tion, and allows for an effective reduction of the degrees of 
freedom of the system. Kadanoff's block-spin rule is prac- 
tically implemented by approximating a 2'^°-spin HE A as 
a 2''°~i-spin HEA. Such an approximation is practically 
performed by imposing that some observables of the 2'^°- 
spin HEA are equal to some corresponding observables of 
the 2'^"^i-spin HEA. For large fco, the method is asymp- 
totically exact, and so its predictions on the critical fea- 



tures of the system. The method has been tested in the 
simplest case of Dyson Hierarchical Model, which is the 
ferromagnetic version of the HEA, and the resulting pre- 
dictions on the critical exponents are in good agreement 
the results obtained heretofore ]2]. 

The method has been then applied to the HEA, and 
identifies the existence of a phase transition in the sys- 
tem, yielding a prediction on the critical exponent ly 
related to the power-law divergence of the correlation 
length at the critical point. Above the upper critical 
dimension a < 2/3, the results for v are in very good 
quantitative agreement with those given by the replica 
method [3] even for small ko = 3, 4. Below the upper- 
critical dimension a > 2/3, the e-expansion for the crit- 
ical exponents performed within the replica method is 
not predictive, since the first few orders have a non- 
convergent behavior, and higher orders are not known. 
Hence, a quantitative comparison between the real space 
approach and the replica method is not possible. On the 
contrary, Monte Carlo (MC) simulations ]10] for a diluted 
version of the HEA yield a prediction for the critical ex- 
ponents in this region. These are in disagreement with 
those of the real space approach. This discrepancy could 
be due both to the smallness of ko in the real space ap- 
proach, or to the non-universality of the exponent v when 
passing from the model on the hierarchical lattice to the 
diluted model, or to the fact that correlations between 
the spin couplings have been neglected in the real space 
approach. Accordingly, the quantitative estimate of v 
below the upper critical dimension is a still untamed is- 
sue, which could be suitable for future investigations and 
developments of the present real space method. 
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